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We propose a method for Monte Carlo simulations of systems with a complex 
action. The method has the advantages of being in principle applicable to any such 
system and provides a solution to the overlap problem. We apply it in random 
matrix theory of finite density QCD where we compare with analytic results. In 
1 this model we find non-commutativity of the limits /i — > and N — > oo which 

could be of relevance in QCD at finite density. 



1. INTRODUCTION 

There exist many interesting systems in high energy physics whose action 
contains an imaginary part, such as QCD at finite baryon density, Chern- 
Simons theories, systems with topological terms (like the (9-term in QCD) 
and systems with chiral fermions. This imposes a severe technical problem 
in the simulations, requiring an exponentially large amount of data for 
statistically significant measurements as the system size is increased or the 
critical point is approached. Furthermore, the overlap problem appears 
when standard reweighting techniques are applied in such systems and it 
becomes exponentially hard with system size to visit the relevant part of 
the configuration space. In Ref. [1] it was proposed to take advantage of a 
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factorization property of the distribution functions of the observables one 
is interested to measure. This approach can in principle be applied to any 
system and it eliminates the overlap problem completely. In some cases it is 
possible to use finite size scaling to extrapolate successfully to large system 
sizes where it would have been impossible to measure oscillating factors 
directly. The method has been applied successfully in matrix models of non 
perturbative string theory (IKKT) 1 , random matrix theory of finite density 
QCD (RMT) 2 as well as the 2d CP 3 model, the Id antiferromagnetic model 
with imaginary B and the 2d compact U(l) with topological charge 3 . In 
this paper we present our results for RMT which we study in order to test 
the factorization method against known analytical results. We also discuss 
an observed non-commutativity in the limits p — > and iV — > oo which 
maybe relevant to Taylor expansion and imaginary p approaches to the 
problem of finite density QCD. 

The factorization method has the important property that it can be 
applied to any system featuring the complex action problem. Let a system 
be given by a partition function Z = J dAc~ So e tT and the corresponding 
phase quenched model Zq = J dAe~ s ° where S = So — iT is the action of 
the system with its real and imaginary parts. A represents collectively the 
degrees of freedom of the model and in our case it corresponds to a set of 
N x N matrices. In case we are interested in measuring some observable 
O, we consider the distribution functions po( x ) = — O)) and p^\x) = 
(S(x — O))o, where (. . .) refers to Zq. Then we define the fiducial system 
Zo,x = J dAc~ So 5(x — O), the weight factor wq{x) — {e lT )o,x and the dis- 
tribution po{x) factorizes po{ x ) = c Po( x ) w o( x ) where C = (c 4r ) . Then 
(O) = dx x p^j (x) wq (x) . The 5-function constraint is implemented 

in our simulations by considering the system Zo,v — J dAe~ So e v (°^ where 
V(z) = \"l{z — £) 2 and 7, £ are parameters which control the constraining 
of the simulation. The results are insensitive to the choice of 7 as long as 
it is large enough. Then we have that wq(x = (O)j.y) = (c tT )iy. The 
distribution of O in Ziy has a peak x and the quantity V'(x) is the value 
of f@\x) = jMnp^(x) at x = x. The function Pq\x) can be obtained 
by integrating an analytic function to which we fit the /q \x) data points. 

By applying this method we force the system to sample configurations 
which give the essential contributions to (O), something that would be 
exponentially difficult with system size in the phase quenched model, elim- 
inating this way the overlap problem. This already allows us get close to 
the thermodynamic limit with modest computer resources. Furthermore we 
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obtain direct knowledge of Wo(x) and po{%) which allows us to understand 
the effect of Y. This is important for understanding the properties of the 
system when T plays a crucial role. 

2. RMT OF FINITE DENSITY QCD 

We consider RMT with one quark flavour and zero quark mass 4 . The 
model is chosen in order to study the correctness and effectiveness of the 
factorization method, since one can compare results with known analytical 
solutions even for finite N. The observable we measure is the "quark num- 
ber density" v as a function of the chemical potential p, and we consider 
the distribution functions pi (x) , where i = R, I corresponds to the real and 
imaginary parts of v respectively. Notice that the effect of T is dramatic, 
causing a discontinuous transition in v. Our results 2 nicely reproduce the 
exact results known for finite N and we are able to achieve large enough 
values of N < 48 to obtain the thermodynamic limit. In Figure 1 we show 
the plots of the distribution functions prj(x) for p = 0.2. Unfortunately, 
the function wr(x) is not positive definite and the important contributions 
come from the region where it changes sign. As expected, we find that 
finite size scaling does not work as well as in the case of the IKKT model 
1 (although we obtain agreement up to order of magnitude for the values 
of N < 96 that we explored). We also find it very difficult to explore the 
critical region near the phase transition point p c = 0.527. . . for N > 8 
since |iUj(x)| becomes very small. Since RMT is a schematic model of finite 
density QCD, we expect that the factorization method will be useful to 
explore the phase diagram of QCD. 

In our simulations we find that for certain observables the limits p — > 0, 
N — > oo and N — > oo, p — > are not equivalent. In real QCD, the former is 
easy but not the latter. That this situation is possible can already be seen 
at the partition function level where Z(p, N) = c K [l + - — ^ — j(N + 1, n)], 
k = —Np 2 is equal to 1 and respectively, but it turns out that the same 
is true also for observables like -§-{vr)q and wr(x). (vr), however, is well 
defined in this limit as expected. In the thermodynamic limit (vr)o = p 
for < p < 1. In Figure 2 we see that this limit is approached like 
(vr)o — p ~ 0(l/N) but only if p > p c (N). We find that the value of p c 
is consistent with p 2 . <~ 1/N. A circle with radius p c — l/y/N contains 
only one eigenvalue of the matrix on average. For N -C we nn d that 
(vr)o = 0. This can be seen clearly from the second plot of Figure 2, 
where the distributions Pr\x) for N — 8 peak around zero for p < p c 
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and their peaks get closer to jj, as fj, becomes larger than fj, c . Therefore 
limjv^oo lim M _>o ^( u r)o = ^ Hm^o lini/v^oo ^(vr)q = 1. Similarly 
we find that limjv^oo lim^^o wr{x) = ^ lim^o limjv->oo wr(x) = 1. 
Details will be reported elsewhere. 
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Figure 2. The difference of {up)o from its A 1 " = oo value and p^' (x) for N 
various values of fi. 
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